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This work addresses the resurgent properties of the cusp anomalous dimension’s strong coupling 
expansion, obtained from the integral Beisert-Eden-Staudacher (BES) equation. This expansion is 
factorially divergent, and its first non-perturbative corrections are related to the mass gap of the 
0(6) cr-model. The factorial divergence can also be analysed from a resurgence perspective. Build¬ 
ing on the work of Basso and Korchemsky, a transseries ansatz for the cusp anomalous dimension is 
proposed and the corresponding expected large-order behaviour studied. One finds non-perturbative 
phenomena in both the positive and negative real coupling directions, which need to be included 
to address the analyticity conditions coming from the BES equation. After checking the resur¬ 
gence structure of the proposed transseries, it is shown that it naturally leads to an unambiguous 
resummation procedure, furthermore allowing for a strong/weak coupling interpolation. 


I. INTRODUCTION AND SET-UP 

The cusp anomalous dimension plays a central role in 
the study on many observables in four dimensional gauge 
theories. In supersymmetric Af = 4 Yang-Mills theory 
(SYM), it appears when studying the scaling behaviour 
of the anomalous dimension of a Wilson loop with a light¬ 
like cusp in the integration contour, in the SL(2) sector 
of the theory [1[. The Wilson loop operators carry a 
Lorentz spin S and a twist L, and for the case of large 
spin and L ~ In S the scaling behaviour of the minimal 
anomalous dimension is [2-5] 

7 s,l ( 9 ) = (2T cusp ( g) + e(g,j))\nS + O (L -1 ) , 

where (47rg) 2 = A = <?ym-^ is the’t Hooft coupling. Also, 
j = L/lnS is the only dependence on the twist from 
the leading contribution to the scaling, and taking j = 
0 leaves us with the twist independent cusp anomalous 
dimension [6, 7]. It is a function solely of the coupling, 
and has been thoroughly studied in different regimes. 

At weak coupling this function can be expanded in 
powers of g 2 , with coefficients determined from perturba¬ 
tion theory [8, 9], and the corresponding series is conver¬ 
gent. At strong coupling, through the AdS/CFT corre¬ 
spondence [10], one can obtain an expansion in g~ 1 from 
the semiclassical analysis of the energy of folded spinning 
strings in AdS.5 x S 5 , where the Lorentz spin and twist 
become angular momenta of the string solution [11, 12]. 

Studying the interpolating region between weak and 
strong coupling is difficult, and integrability played a cru¬ 
cial role. The all-loop Bethe ansatz for J\f = 4 SYM [13- 
15] led to a set of integral equations, the BES equations 
[4, 16-18], describing the anomalous dimension, valid for 
any arbitrary coupling (with the FRS equations [4] valid 
for any scaling parameter j). In terms of an auxiliary 
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function 7 (2 gt), the BES equation can be written as 

= K (2 gt, 0) - 2 g [ -^—K(2gt, 2gt') 1 (2gt'). 

2 gt Jo e* — 1 

where the BES kernel K (t, t') can be found in [16, 17, 19]. 
This auxiliary function is related to the cusp anomalous 
dimension by 


Tcusp { 9 ) 


= 8 lim 

t—y 0 


l(2gt) 
2 gt 


Solving these equations at weak coupling returned 
higher terms of the convergent perturbative expansion 
for the small j <C 1 region [17]. For intermediate cou¬ 
pling g ~ 1 a smooth solution to the BES equation was 
found numerically [20]. At strong coupling different at¬ 
tempts were made at solving the BES equations [21-23], 
and in [24, 25] a solution was found leading to a strong 
coupling expansion. This approach consisted in noticing 
that a change of variables 


ft (t) = 7 (t) ^1 + cosh , 

returns a simpler set of coupled integral equations for 
fl(t), which are then solved using Fourier methods [26, 
27]. One subsequently obtains a solution of the BES 
equation in the form [19, 28] 


ft (it) = fo (t) V 0 (t) + fi (t) Vi ( t ), 
where (x = 8irg) 

.fo ( t ) = —1 + 

+ X 2t 


(1) 


n> 1 


, ,urm 

C + [n ’ X) + c - <"• x) 


h m s Y. 


TlX 


T> 1 


TT+ (a*} tj- ( 

c+ {n,x) U ^M^+C-(n,x) U ° ( 2 } 


nx — 2 1 


nx + 2t 


The functions Uq 1 and Vqj can be written in terms of 
Whittaker functions of 1 st and 2 nd kinds, but for our pur¬ 
poses we only need their asymptotic expansions for large 
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x , which can be found in Appendix A. The coefficients 
c±(n,x) are determined from analyticity conditions on 
the solution f l (it) (given that 7 (it) is an entire func¬ 
tion): from the expresions for f n (t) it already has the 
correct pole structure, but one still needs to impose the 
existence of zeroes at 

^zeroes = = | (t - ^ , i G Z. (2) 


This condition can be re-written as 


i = E 


C+ (n,x) 
n- ae 


n + a? 


%(f)c« + [ y(f )», W ) 

(r r (f) O( + [/ 0 -(f )»r (a,)),! 


Two important questions still remain at this point: are 
we aware of all of the non-perturbative phenomena defin¬ 
ing the analytic properties of the cusp anomalous dimen¬ 
sion? How can we systematically deal with a non-Borel 
summable asymptotic series? To answer both these ques¬ 
tions we will now turn to the theory of resurgence. 

Resurgent functions have been seen in a wide range of 
systems. In mathematics they appear for example as so¬ 
lutions of differential and finite difference equations (see 
e.g. the well studied cases of Painleve I, II and Riccati 
non-linear differential equations [31-34]). Analogously, 
often one can only determine physical observables in spe¬ 
cific regimes of the coupling of the theory via a series 
expansion such as 

{0(g))^Y,°^~ k - ( 5 ) 

fc >0 


where r (a) is the ratio of functions V\(nx/2) and 
Vq (nx/ 2), as defined in Appendix A. This analyticity 
condition allows us to determine the coefficients c± (n, x) 
order by order as expansions in large coupling x. Once 
these coefficients are known, the cusp anomalous dimen¬ 
sion is given by (x = 8ng) 


However, these expansions are often asymptotic: the co¬ 
efficients are factorially divergent, with large order be¬ 
haviour 


r (k + p) 

A k +P ’ 


k > 1 . 


( 6 ) 


= l-2/i(0) (4) 

= !- 2 E [ c +( n > x ) u o (^r) + c_(n, x)Uq (^pj . 

n> 1 

The strong coupling expansion found in this way is 
asymptotic. Moreover, the series is non-Borel summable 
for positive real coupling, due to the existence of singu¬ 
larities on the positive real axis of the Borel plane, which 
give rise to non-perturbative, exponentially suppressed 
corrections at strong coupling. In order to understand 
the analytic properties of the solution to the BES equa¬ 
tion at strong coupling, one needs to account for all the 
non-perturbative phenomena in this limit. In [19, 28] the 
above procedure was taken a step further and the per¬ 
turbative coefficients around the first non-perturbative 
correction were determined. 

Both scaling function e(g,j) and cusp anomalous di¬ 
mension T cusp ( g ) have non-perturbative corrections. In 
[3] it was proposed that the scaling function e(g,j) at 
strong coupling is directly related to the energy density of 
the ground state of the 0 ( 6 ) non-linear er-model embed¬ 
ded in AdSs x S 5 (taking j /2 to be the particle density). 
Moreover, the non-perturbative corrections appearing in 
e (g,j) at strong coupling are given by the mass scale 
(mass gap) of the 0(6) model. Agreement between these 
two quantities was checked in [19, 29, 30], at the level of 
the first non-perturbative correction to the scaling func¬ 
tion. As for the cusp anomalous dimension, as it solves 
a different integral equation altogether, such a relation 
was less expected. Nevertheless, in [28], it was shown 
that the first non-perturbative correction to the anoma¬ 
lous dimension is exactly given by the square of the 0 ( 6 ) 
mass gap. 


A, (3 are numbers related to the position and type of sin¬ 
gularities of the related Borel transform. 

It is well known that this divergence hints to the ex¬ 
istence of non-perturbative phenomena unaccounted for 
in the perturbative series expansions. In physical set¬ 
tings, the existence of non-perturbative phenomena has 
been long noticed in the contexts of quantum mechan¬ 
ics [35, 36] and quantum field theories [37], associated 
to instantons [38] and renormalons [39]. In these exam¬ 
ples, the existence of asymptotic multi-instanton sectors 
allowed for a complete unambiguous description of the 
energy eigenvalues via a transseries solution and resur¬ 
gence [40-42]. Since then, the asymptotic behaviour of 
perturbation theory and the resurgence behind it was 
seen to exist in many different examples in physical sys¬ 
tems, from quantum mechanics [43-45], to large N gauge 
theories [33, 34, 46-53], quantum field theories [43, 54- 
60] and topological strings [61-63]. 

To account for all non-perturbative phenomena, one 
upgrades our perturbative expansions into a transseries 
[50]: a formal expansion in both perturbative variable g 
and non-perturbative monomials e~ Ag . Schematically 

(O (g, a)) = £ a"e-^V" E OfcV*, (?) 

n>0 k >0 

where a is a parameter to be fixed from some boundary 
conditions specific to each problem. The transseries is 
a formal object, since for each non-perturbative sector 
labeled by n one has an associated asymptotic expan¬ 
sion 4>„ {g) ~ g@ n X]fc>o with coefficients grow¬ 

ing factorially at large orders. However, these sectors are 
not independent of each other: they are sectors of a resur¬ 
gent transseries, whose large-order growth is intimately 
related [33]. A resurgent transseries is an expansion like 
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(7), where the coefficients of one sector O^'" 1 are related 
to, i. e. resurge in, the coefficients of neighbouring sectors 
C>| m 1 (to close to in'). For example, for the perturbative 
sector (n = 0 ), one expects a direct relation to the n = 1 
sector 

°k ] ~ F %+l n) + fa* + •••), k » 1. (8) 

The exact expressions for these large-order relations [33] 
can be determined via resurgent analysis (for a introduc¬ 
tion to resurgence see [33, 64-66] and references therein). 
The associated Borel transforms B [$„] (s ), 1 have a non¬ 
zero radius of convergence and singularities on the corre¬ 
sponding Borel s-plane at positions s = nA, n £ N. 

At this point we have a formal solution for our observ¬ 
able, and we now need to retrieve physical information 
from the asymptotic series <f>„(g). This is done via Borel 
resummation: the calculated Borel transform has a non¬ 
zero radius of convergence, and one can determine the an¬ 
alytic function associated with each series B [$„] (s), ei¬ 
ther exactly or by finding an approximate analytic result 
via the so-called Borel-Pade approximants [33, 50, 67]. 
Once the function or its approximant is known we then 
perform a resummation via a Lapace transform 

r+oo 

S<t> n (g) = / dse- S 9 S[$„](s), (9) 

Jo 

and the full answer for the observable is given by the 
transseries with each of its sectors resummed. This can 
only be performed if no poles exist in the direction of 
integration on the Borel plane, in this case on the positive 
real line. If instead A is positive and real, the positive real 
line is called a Stokes line (singular direction on the Borel 
plane), and the series is said to be non-Borel summable: 
only lateral resummations can be defined: 

z’+oo e^ le 

S±$n(g)= dse- ss Bi<S> n ](s), (10) 

Jo 

these lateral resummations differ by a non-perturbative 
ambiguity (5+ — <S_) $„ (g) ~ e~ Ag , which is purely 
imaginary when the coefficients O'^ 1 are real and the 
Stokes line is along the real axis. Now the importance 
of having a resurgent transseries becomes apparent: due 
to the relations between different sectors, by taking into 
account the full transseries and a specific value for the 
transseries parameter a, the ambiguities between differ¬ 
ent sectors cancel each other, and one is left with a non- 
ambiguous real-valued result. This is called median re¬ 
summation (see [ 68 ] and references therein). 

Recalling the transseries (7), let us assume that the 
positive real axis is a Stokes line and choose the lateral 


Borel resummation <S + for every sector of the transseries. 
This resummation will have real and imaginary parts 

5+ < I> n = — (5+ + iS_) 4>„ + - (5+ — S-) 

= ‘Sft < f> ri + idz/d’n, (11) 

The imaginary contribution Si& n is just the ambiguity 
coming from the sector 4> n . We can now determine the 
real and imaginary parts of the resummed transseries [ 68 ] 

S + (0(g,a)) ee 

n>0 n>0 

= S fl < 0 > + iSj< 0 >. 

The ambiguity in the resummation of the transseries is 
just its imaginary part (<j = ctr + i oy) 

5/ (O) = Im Re + a R Im (V 1} ) + 

+2a R (JiRe ( F (2) ) + (a 2 R - aj) Im (f< 2 )) + • • • . (12) 

Median resummation is a specific prescription to cancel 
this imaginary contribution to the resummation of the 
transseries along the positive real axis to all orders, by 
some carefully chosen values of the transseries parameter 
a = do - 2 This cancelation happens to all orders, and we 
are left with a real unambiguous answer 

Sr (0(g, cr 0 )) = Re (V 0 )) + do^Re (F 1 ) + (13) 

+ (ofU - °o,/) Im ( f(2) ) +■■■ ■ 

This resummed result can then be interpolated from the 
original regime where the asymptotic series were defined, 
into any complex value of the coupling g, taking into 
consideration any crossing of singular Stokes lines. The 
systematic resummation and interpolation from asymp¬ 
totic series using resurgence has recently been addressed 
for different problems [53, 59, 63, 69, 70]. 

The aim of this paper is to perform a resurgent anal¬ 
ysis of the expansions found in [28] for the strong cou¬ 
pling regime of the cusp anomalous dimension. We start 
from the solution to the BES equations presented above, 
and enforce the analytic properties at the level of the ex¬ 
pansion coefficients. We then determine the structure of 
singularities of the Borel transform associated to the per¬ 
turbative sector by means of a Borel-Pade approximant. 
This allows us to finally propose a transseries ansatz for 
the cusp anomalous dimension which encompasses all the 
expected non-perturbative phenomena existing at strong 
coupling. Using this ansatz in the analyticity conditions, 
we determine the coefficients of our transseries, solved 
order by order for every sector. 


1 Borel transforms are determined by inverse Laplace transforms 2 In simple cases it was seen that icrj = -4 was enough to cancel 
to each term in the expansion, or equivalently g~ k —t s fc- 1 /r(fc). the ambiguity [68], with residual freedom left in the real part or. 
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Equipped with the series expansion for perturbative 
and non-perturbative sectors of the cusp anomalous di¬ 
mension, we then check its resurgent properties via the 
large-order relations. Along the way we determine the 
relevant Stokes constant associated with the Stokes tran¬ 
sition across the positive real line. 

We finish by using the methods of median resumma¬ 
tion to systematically define a non-ambiguous resummed 
result valid at any value of the coupling, which encodes 
the analytic properties of the solution and can be used to 
interpolate between strong and weak coupling regimes. 


II. SINGULARITY STRUCTURE OF THE CUSP 


In the interest of finding the correct transseries solu¬ 
tion for the cusp anomalous dimension, we first anal¬ 
yse its perturbative asymptotic series. With that goal in 
mind we assume the coefficients c± (n, x) have a simple 
(asymptotic) expansion in powers of x where x = 8ng: 


+ OO 

c± (n,x) = .W 0 )^ 1 / 4 ^^) (n) x~ k . (14) 

k =o 


Substituting this ansatz into the analyticity conditions 
(3), and making use of the asymptotic expansions in Ap¬ 
pendix A, as well as properties of sums found in Ap¬ 
pendix B, one obtains for each power in x~ k relations 
for the coefficients </4°’ ± ' , (n) depending on the ones for 
lower k. Solving these relations iteratively (in the same 
way as was done in [28] for the first few coefficients), we 
determine 


^• ±, (») = E «£? 

m—0 


-nAO,±) 


n' (+> ( 15 ) 


where the expressions for </>Q 0,± ' ) (n) are very simple and 
can be found in Appendix B. The analyticity condition 
imposes restrictions on the coefficients /3±(0): 


Afc(O) = ±1/4. (16) 


The general solution for the numerical coefficients 
is then simply given by 


4 Q 


(o,+) 

k,k 


4 Q 


( 0 ,-) 

k,k 


k -1 


k—r 


-£g& +) £ m K. 


( 0 ,+) 


r =0 
k -1 


£=0 

k—r 




r —0 


£=0 


(17) 


for to = k, and for 0 < a < k 

2<?i°; +> = -4 E 

m=a -\-1 

a r k 


(18) 




r—1 £=0 
k k 


- E EoKl-Ar.’E"^’ 


r—a+1 m—r 
k k 


£=0 


/ j / j ^k—r,m—r lv m—a lvl r— l,r—a— 1 ’ 


r—a+1 m—r 


where similar solutions can be written for Q^l’ ■* by 
exchanging Km +) -H> Km~\ and Af£°’ +) 

The definitions for the coefficients 

m(™M ail( ] can be found in the Appendices. One 

can now determine several coefficients in the expansions 
(14), which was done numerically up to k = 200. For the 
cusp anomalous dimension one then uses the expansion 
(4), and the expansions for the functions Uj~ (^) , k = 
0,1 given in Appendix A. 

The strong coupling perturbative expansion for the 
cusp anomalous dimension becomes 

r cusp (ff) = 2 g (l + T<°> (4 ttp) +0 (e-* 4 "®)) , (19) 

where the perturbative asymptotic expansion is 

+ 00 

r (0 +) = — 2 / liPert (o) ~ ^ (|) r< 0) . ( 20 ) 

k =1 

The expansion coefficients are simply given by 

k —1 m 

r (°) „ 2 1 - fc ^^5 + ( S )Qi°:tl s>m _ s ^+(21) 

m—0 s—0 
k— 1 m 

+ 2l_fe EE 5 ' (s)Qk-~l-s,m- S K™+l\ 

m—0 s—0 


with 


S+(s) 


S-(s) 


r(5/4 + s)T(l/4 + s) 

T (5/4) T (1/4) s! 1 1 J 

r (3/4 +a) r (-1/4 + 3) 
r (3/4) T (-1/4) a! 1 


One can now see that the coefficients grow facto- 
rially for large order k. In fact they grow as T(k — 1/2), 
in agreement with the factorial growth found in [24]. As¬ 
sociated with this factorial growth there will be non- 
perturbative phenomena dictating the asymptotic na¬ 
ture of the series. This non-perturbative phenomena is 
most naturally represented as singularities on the com¬ 
plex Borel plane: we expect to find singularities at posi¬ 
tions s = n A, where A = 1/2 - these will be associated to 
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non-perturbative exponentially suppressed contributions 
of the type e~ nA 47r A 

To study the singularity structure associated with the 
perturbative expansion we determine its Borel 

transform via the usual approach 


B 


p(o) 


(s) = 


+ OO 


p(°) q k 
1 fc +1 * 


^o r(fc + 1) ' 


(23) 


This expansion will now have a non-zero radius of conver¬ 
gence and we approximate the corresponding function via 
the method of Pade approximants: using a diagonal ap- 
proximant or order N = 100 (half the order of coefficients 
calculated for the original series), we determine the best 
fit for a ratio of two polynomials BP^ (r^ 0 )) and analyse 
the structure of poles for this function. This allows us to 
see the position and type of nearest singularities to the 
origin on the Borel plane: in particular, a condensation 
of poles hints to the existence of a branch cut. In Figure 
1 the structure of poles of the Borel-Pade approximant 
is given. We find the expected pole ats = A = l/2, but 
we also find another type of singularity on the negative 
real axis at s = —4 A = —2. 

The first type of singularity had already been known, 
it is directly related to the square of the mass gap of the 
0(6) cr-model embedded in AdSs x S 5 [28]. The fact that 
it lies on the positive real axis prevents us from defining 
a resummation on this axis: we can only define lateral 
resummations (10) which differ by an imaginary ambigu¬ 
ous contribution. However, the second type of singular¬ 
ity lies on the negative real axis. Even though it will not 
give any ambiguous contribution to a resummation for 
real and positive coupling g, the analyticity conditions 
will not be blind to it. These types of singularities have 
been found before in the study of the Painleve I and II 
equations [32-34]. 

If the perturbative expansion for some observable is 
asymptotic, one should upgrade the solution to include 
non-perturbative sectors, into what is called a transseries. 
The most important ingredient to writing a transseries 
solution fully describing our observable is to include all 
possible sectors associated with singularities on the Borel 
plane. In the present case, this means upgrading the 
expansions for c± (n, g ) to a transseries including both 
types of singularities found. 


III. TRANSSERIES AND ANALYTICITY 
CONDITIONS 


We can now write a trasseries ansatz for our coefficients 
c± (n, g) appearing in (4). 

+oo 

c± (n, X] (Tj) = ]T <Vr 2 e (4m2_mi)Af 4 mi V,z), 

mi ,7712 —0 

(24) 

where cr, are the transseries (or instanton counting) pa¬ 
rameters, A = I and c± l \n, x) are perturbative ex- 


PolesofBP™!™ 


Figure 1: Poles of the diagonal Borel-Pade approximant of 
order 100 for the perturbative series of r(2<?)/(2<?) — 1. There 
accumulation of poles in both positive and negative real di¬ 
rections, starting at s = A = 1/2 and at s = —4 A = —2. 
Note the existence of spurious poles away from the real line: 
non-stable numerical effects of the Pade method, which move 
away by choosing different non-diagonal approximants. 


pansions around each non-perturbative sector ( 7711 , 77 x 2 ): 
c±\n,x) are simply the perturbative expansions (14) 
found in the previous Section, while the others will gener- 
ically be 


+OO 

c£ mi) (n, x) ~ x^ mi) Y, ^" fc 4 m ’ ±) (n) n ±1/A , (25) 

fc =0 

where /3± (rn) are numerical factors associated with the 
type of branch cuts on the Borel plane. In general one is 
interested in the transseries solution at some particular 
value of the coupling x (usually real positive), and the 
parameters <Ji are fixed by some physical input. More¬ 
over, having a non-zero cr 2 for positive real x would lead 
to unstable, exponentially enhanced contributions, which 
physically we know should not be included. Nevertheless, 
such sectors need to be introduced in order to account for 
all the analytic properties of the problem at hand, even if 
in the end, once resummation is performed, the parame¬ 
ter cr 2 gets fixed to zero. 

Another important issue is that of resonance [32— 
34, 66 ]: for transseries with more than one type of ex¬ 
ponential behaviour e -"‘di- m 2 ^ 2 ] if there are values 7777 
such that m\Ai + 777 2 A 2 = 0 , a phenomena called res¬ 
onance occurs. In the case of transseries solutions to 
non-linear differential equations, it is seen that the asso¬ 
ciated recursion relations break down at these locations 
unless we enhance the perturbative expansions (25) 
to include other non-perturbative sectors, such as some 
finite number of powers of log(a:). In the present case we 
will find resonance when 777.1 = 4777 2 , and one can expect 
a rich structure like the one found for the Painleve solu¬ 
tions [33, 34], For our present case we will limit the study 
to the transseries contributions up to "two instantons", 
777 = 2 in (24), and will not reach these structures in our 
analysis. 

Having written the transseries ansatz, we now need to 
substitute it in the analyticity conditions (3). After some 
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algebra, we can re-write these conditions as 


where 


2 = 


+ OO 

E 

rri 2 ,m =0 

LapbJ 

x E 

r—0 


a™a™ 2 e - A %( m - 4m2) x 


(26) 


r —4|a|r ^ ^ + p(.m,r) ^ ^ 


which need to be obeyed for every zero a = £ — ^, £ € Z. 
Throughout this paper, we focus on the contributions 
with to < 2 and to -2 = 0. Therefore we will leave the issue 
of resonance and higher non-perturbative corrections for 
subsequent work. We briefly note that the last sum in 

which will be 


4|a| 




(26) goes up to the integer part 

non-zero if m > 4 | a |. For a# < 0 this already happens 
for to = 1 and ag = —1/4. For ag > aq = 3/4 > 0, 
the first non-perturbative contribution will be at m = 3. 
The case to = 4 is somewhat special, as both ao = —1/4 
and = 3/4 return a non-zero sum (for to > 5 there 
will be two negative zeroes). It would also be at this 
point (taking m 2 = 1) that we find the first instance of 
resonance: it is likely that the two effects will mix in the 
analyticity conditions. 

Taking m 2 = 0, then the functions F± are 3 

F (m,r) (n,x) = T/3+( m - 4 l«k)- 1 /4 


x 1 

+ OO 

*E* 

k—0 

A 


- fc E 

n> 1 


X 

a—4|a|r,+) 


(27) 


(n) 


n — a 


xi 0 ,i (x, -n) 5 r ,0 


1,1 


(x, —n) 


A 0 ,i (x, —n) 


x 1 R r 


(x,a) 


F { ™' r) {n,x) = ^(™- 4 l«k)-3/4 x 

+ OO 

x E 


(28) 


+°° ,(m-4|a|r,-) , , 

-fc LA _W 

fc=0 n> 1 


(n + a) n 

A),o (x, ~n) f c^ E 1 ' 0 + n R r (x, a)] . 

\ A),o [x, —n) J 


We can now solve (26) for each to, and this was done 
for to. = 1,2, in the same manner as for the perturbative 
series in the previous Section. We found that 

/Ml) = ^/Ml) = ^/3i(2) = ±^ (29) 

and also 

4 m,+) («) = 

e=o 

= Eflir/t-’w. (so) 

^=0 



= -P (m) 4°’ +) («), 



= 0, 

(31) 


= P (m) 4 0r) (n-1). 



The numerical coefficients V are also determined: 

(1) ole 3 -/ 4 T (3/4) 

2 T(5/4) ’ 

■p(2) = _^(!)j 2 . (32) 


It is worth noting that these coefficients seem to depend 
on the transseries parameter, which it is not yet fixed. 
The value of the transseries parameters can vary with 
the value of the coupling x: even if we fix them to a 
particular value, when we move on the x-plane, these pa¬ 
rameters will jump in value when crossing any Stokes line 
(lines where there are singularities on the Borel plane) - 
such jump will be governed by the so-called Stokes auto¬ 
morphism. So how can we interpret the numerical values 
'p(m) they depend on a parameter which will change 
its value? In fact, the analyticity conditions are solved 
in a specific direction on the x-plane: the one which was 
chosen to perform the asymptotic expansions. In this re¬ 
gion we will have a specific value <ri = off, and the r p( rn '> 
will be fixed at that position on the x-plane. 

Once the c±‘\n, x) are determined for m = 0,1, 2, we 
can write down the corresponding transseries for the cusp 
anomalous dimension (4). One can write down the full 
transseries solution corresponding to (24), but for the 
present we will take m 2 = 0. Then 


+OO 


r/usp (<?, < j ) — 2$ 1 + / ' 


a m e -mA f p(n 


\x) , (33) 


m =0 


where a = was the perturbative contri¬ 

bution previously calculated (21). The two first non- 
perturbative corrections can be written as (to. = 1,2) 


r( m )(a:) 

where /? (to) = — ™ 

r 


-|-oo 


-k 


+00 

E 

k =0 

and the coefficients are 4 


(m) 


= -4 


■p(n 


(Pd)) ? 


(34) 


(35) 


See Appendices for the expansions used, and Rk ( x , a) 
Rk (x , |a|) if a > 0 while Rk ( x , a) = Rk ( x , |o;|) if a < 0. 


The first coefficients around the m = 1 sector are in agreement 
with [28]. 
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r (m) _ 

1 7. — 


. , -p (m) 
_oi-fc__ v 

" / _ / - s \ m. 


(■pd)) 

EE 5 -( s )^-J-,E 4a+1_r ^ 0 ' 


(36) 


A; a 


\a —0 s—0 
/c— 1 a 


r—0 


+ E E 5 + ( S ) , k > 1 . 


a—0 s—0 


We have now calculated the perturbative coefficients 
around the first two non-perturbative sectors. If we anal¬ 
yse the growth of these coefficients, we find the same fac¬ 
torial growth r(fc- 1/2) for m = 0,1, and the factorial 
growth r(fc + l/2) for m = 2: not only the perturba¬ 
tive series is asymptotic, but so are the non-perturbative 
ones. Moreover the singularities which lead to these sec¬ 
tors lie on the positive real axis, and thus one cannot 
properly define a single integration contour on which to 
perform the resummation of the Borel transform. 

Nevertheless, many lessons have been learned by now 
on the cases of resurgent transseries: if our transseries is 
resurgent, then there is a way of defining a single non- 
ambiguous result which properly cancels the imaginary 
ambiguity at all non-perturbative orders. But in order to 
use these results from resurgence, we first need to check 
that our transseries is indeed resurgent. 


IV. RESURGENCE AND LARGE-ORDER 

Let us now check that the transseries formed by the 
asymptotic series ]?( m ) is indeed resurgent, i.e. that the 
coefficients rj, m \l^ m ' of neighbouring sectors m,m' are 
related. To perform this check we use the coefficients 
r/ m) of the asymptotic series (34) and check if their large 
order behaviour coincides with the large-order relations 
predicted by resurgence techniques. These are relations 
between the large order k 1 r[ m ' ) of one sector with 

the low order coefficients of a nearby sector / 

Take the transseries (33); assuming this transseries is 
resurgent, we can use the so-called alien calculus to de¬ 
termine the discontinuity of each asymptotic series 
across singular directions (or Stokes lines). In our case 
(given we have taken m 2 = 0) we only have one singular 
direction: the positive real axis 0 = 0. Resurgence then 
tells us [33] that the discontinuity of the perturbative 
series along this direction is (recall that A = 1/2) 

Discor*- 0 -* (x) = -s ie ~ A %r (1) - (S /) 2 e ~ 2A %rM + ■■■ . 

(37) 

There is one unknown constant in the above relation, the 
Stokes constant Si. As we will see, the large order rela¬ 
tions will allow us to determine this constant with great 
accuracy. This step is extremely important as the Stokes 
constant plays a crucial role in the ambiguity cancelation 
and resummation. 


From the discontinuity, we can use Cauchy’s theorem 
to determine large-order relations [38]. Schematically, 
one writes 


r (0 >(«) = (f 

J c 0 


UJ — Z 

+00 


dto rW(cu) 

2M w — z 
dui Discor(°)(w) 
27ri u> — z 


(38) 


+ 




In certain conditions, it can be shown by scaling argu¬ 
ments that the integral at infinity does not contribute 
[36, 71]. Expanding the r.h.s for large z, using the resur¬ 
gence relation for the discontinuity, and finally comparing 
equal powers of z for the expansions in both sides of the 
equation, we arrive at the relation 


p(o) 

1 k 


SiV {1) p(i) T (k — \ — h) ~ 

27ri ^ ^4 ^~h~h 

h=0 

( Sl yw) 2 ^ 

2m ^ h (2A) k - 1 ~ h 

h =0 v ’ 


(39) 


k » 1. 


This formula states that if resurgence is expected, then 
the large k behaviour of the perturbative series is dic¬ 
tated by the coefficients of the first non-perturbative sec¬ 
tor, and then, more exponentially suppressed (2~ fe ), the 
coefficients of the second non-perturbative sector appear, 
and so on. The proportionality constant is once again the 
Stokes constant S/. Taking the ratio (which removes the 
dependence on the yet unknown Stokes constant) of two 
consecutive coefficients, and assuming k 1, we have a 
series (asymptotic again) 


do) 


+ OO 


-p(°) A 

1 fc+i 




—h 


(40) 


h —0 


where the coefficients Ch can be predicted from the orig¬ 
inal large order relation (39). The first coefficients are 


1 1 r (1) 

c 0 = 1; ci = c 2 = - + ■ ■ ■ • (41) 

We can now check the convergence of rj. 0 ' to the coeffi¬ 
cients Ch by successively removing the previous coefficient 
from the ratio. For example to check the convergence to 
the coefficient C 2 we analyse 


( (r^ - ^ ~ c °) k ~ Cl ) k C2 + ° ^ ( ' 42 ' ) 

In Figure 2 we present the convergence to coefficient 
C 10 ■ In order for this convergence to be correct, all of 
the previous coefficients need to be correct to a very high 
accuracy, since factorial errors propagate rapidly. In this 
figure, the original ratio (in red) is shown, together with 
two related Richardson transforms which speed the con¬ 
vergence of this series in 1/n (see [32, 49]). The error be¬ 
tween the numerically calculated coefficient (via Richard¬ 
son transforms) and the predicted result from the large- 
order formula is of order 10 _ 7 . 




Figure 2: Convergence of the large order ratio of perturbative 
coefficients to the predicted result related to the first non- 
perturbative sector. In red is shown the original ratio, and 
in blue the corresponding Richardson transforms of order 2 
(light blue) and 6 (dark blue). In light green the predicted 
value for the coefficient Cio is shown. 


Figure 3: Convergence of the perturbative coefficients to the 
Stokes constant. In red the original series is shown, and in 
blue the second Richardson transform with increased conver¬ 
gence. In dark green the calculated value for the Stokes con¬ 
stant is shown. 


If instead of dealing with the ratio of coefficients we 
analyse the following 


2 nA k -5 rj: 0) 
r(fc-§)4 1} 


LSipW Y, 


+oo -p(l) 


IT' r (k - i 


■p(i) 
h =0 1 0 


r(fc- 


iSiP™ + O (k~ 



(43) 


we directly obtain a convergence to the unknown Stokes 
constant. In Figure 3 this convergence is shown, with 
both the original series and a related Richardson trans¬ 
form. The increase in convergence speed from the 
Richardson interpolation method allows us to determine 
the Stokes constant to a very high accuracy. Up to an 
error of 10 -23 we find 


iSiPW 


IT (3/4) 
2 T(5/4) ' 


(44) 


The same ideas were repeated for the asymptotic series 
of T 1 ), whose large order will be directly related to co¬ 
efficients of T 2 ), finding that resurgence predictions also 
worked in this case. Note that unlike the previous cases, 
the large-order behaviour of T^ 2 ) is dictated not only by 
T® but also by T^ (the two nearest singularities on the 
Borel plane will be equally distant from the origin, at 
s = ±A). We conclude that indeed the transseries for 
the cusp anomalous dimension is resurgent, and thus we 
can apply the methods of ambiguity cancelation known 
to exist for resurgent transseries. 


V. AMBIGUITY CANCELATION AND 
INTERPOLATION 


write a non-ambiguous result? And if this is possible, 
can we then interpolate from the strong coupling asymp¬ 
totic expansion to small coupling? 

The answer to the first question is simply yes. The fact 
that we have a resurgent transseries directly tells us how 
to obtain a non-ambiguous result even when resumming 
in directions which are non-Borel summable, i.e., along 
Stokes lines. 

In order to verify the ambiguity cancelation it is suffi¬ 
cient to check that the imaginary part of the resummed 
lateral transseries cancels to higher and higher orders. 
The ambiguity coming from the perturbative series, i.e. 
the imaginary contribution from performing a lateral re¬ 
summation, is of order e~ A4 '"' 9 for each value of the cou¬ 
pling g. This is exactly the order at which the first non- 
perturbative sector starts contributing. In fact, if we add 
the two with the proper choice of parameter er in (7) we 
can see that the newfound imaginary part will now be of 
order e~ 2A4 ' rrg - the order of the second non-perturbative 
sector. In order to see this cancelation, and to determine 
the value of the parameter a = <JiP (yl ' > which brings about 
the cancelation, we first need to resum our asymptotic se¬ 
ries. 

The method for resumming our series is the so- 
called Ecalle-Borel-Pade resummation method. It con¬ 
sists in calculating the Borel transforms for each sector 
B [T( m )] (s), then determining a Pade approximant for 
each Borel transform, and subsequently perform a lat¬ 
eral Borel resummation (we have chosen <S + as in (10)) 
for different values of the coupling 0.1 < g < 4 in order to 
obtain a resummed result for the full transseries for gen¬ 
eral values of g. 5 The imaginary part of the transseries 


The next two questions which follow are: if our 
transseries is resurgent, can we use this knowledge to 


5 It is very important to keep track of the first few terms of the 
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Cancelation of Imaginary Part of 5 + ng) 



Figure 4: Order of magnitude 10~“ (where a is shown on the 
y-axis) of the imaginary part of the transseries for the cusp 
anomalous dimension, when one adds only the perturbative 
series (dark brown), if one adds also the first non-perturbative 
correction (light brown) and if one takes all three calculated 
contributions (light purple). 


for the cusp 

^-1 = 5 + r (0) +ae-5 47r9 l S + r (1) +a 2 e- 4,r9 5 + r (2) +- • • 

2 g 

(45) 

is given by (12). In Figure 4 the order 10 a of the imag¬ 
inary part can be seen if we only include 5+F^ (dark 
brown), or if we include both perturbative and first non- 
perturbative contributions (light brown), and finally if 
we include all calculated contributions (up to second or¬ 
der, in purple). The cancelation is shown for a range of 
values of the coupling, ranging from weak to strong cou¬ 
pling. For example for g = 1.5 if we include only the 
perturbative part, the imaginary ambiguity is of order 
10~ 5 . Including the first non-perturbative sector cancels 
this imaginary part to 10~ 10 . Including all three sectors 
cancels the imaginary part up to 16 decimal places. The 
value of the transseries parameter used was 


o-q = - (1 + i) 


LS'ipW 

2 


1 .-3*1/4 r (3/4) 

2V2 r(5/4)’ 


(46) 


with the Stokes constant given by (44). This is in com¬ 
plete agreement with the proportionality constant ap¬ 
pearing in front of the non-perturbative corrections in 
[28]. 

We are now ready to answer the second question posed 
in this Section: once we have managed to write down a 
real unambiguous transseries result, given by (13) with 


series which might not be included in the Borel transform, and 
thus need to be added at a later stage. Also the overall factor 
can be left out until after the resummation. 



Figure 5: Resummation of the transseries result for differ¬ 
ent values of the coupling, shown in red. For large coupling 
the behaviour is dictated by the perturbative series alone, 
as shown by the dashed line - the truncated summation of 
the perturbative expansion F^, which diverges for small g. 
In this regime another line is shown in blue, corresponding 
to the small coupling expansion (up to 7 loops) of the cusp 
anomalous dimension. The red dots interpolate between these 
two regimes, starting to diverge around value g = 0.2. 


the parameter do determined in (46), can we use this re¬ 
sult to interpolate between strong and weak coupling? In 
Figure 5 we show the truncated asymptotic series result 
in a dashed green line: this result is accurate at strong 
coupling, but diverges for weak coupling. In blue we plot¬ 
ted the weak coupling expansion for the cusp anomalous 
dimension as determined (up to 7 loop order) by [17]. 
Naturally this result diverges for large values of the cou¬ 
pling. In red we show the resummed result including 
the perturbative series and the first two non-perturbative 
sectors. We see clearly that the red dots follow both the 
strong and weak regime closely, starting to diverge for 
g < 0.2. In order to obtain more accurate results af¬ 
ter this point, we would need to include the next non- 
perturbative order. 

In summary, once we have established resurgence of 
the transseries, we are able to write down an unambigu¬ 
ous result, and use that same result to reach values of 
the coupling as small as g = 0.2. Moreover, the infor¬ 
mation encoded in the transseries solution for the cusp 
anomalous dimension goes beyond the strong/weak in¬ 
terpolation for positive real values of the coupling. The 
resummation can be performed for g £ C: one can obtain 
a solution for any value of coupling, having in mind that 
to reach certain values one might have to cross a Stokes 
line and the transseries goes through a so-called Stokes 
transition (in a Stokes transition the transseries parame¬ 
ters will jump in value, and this jump is dictated by the 
Stokes constants and can be calculated from resurgence 
techniques). In other words, the transseries solution en¬ 
codes the analytic properties of the cusp anomalous di¬ 
mension as a function of the coupling. 

Another extremely insightful example of how the 
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transseries encoded the analytic properties of the ob¬ 
servable was recently achieved in the context of matrix 
models [53], where from a large N asymptotic expansion, 
resurgence and resummation techniques allowed the au¬ 
thors to reach not only finite values of the rank N, but 
they were also able to analytically continue their results 
to any complex value of N, once Stokes phenomena were 
taken into account. 


VI. CONCLUSIONS/OUTLOOK 

In this work we presented a thorough analysis of the 
resurgent properties of the cusp anomalous dimension’s 
strong coupling expansion. When analysing the pertur¬ 
bative series we found that there are two types of singu¬ 
larities on the Borel plane in both the positive and nega¬ 
tive real axis. These need to be taken into account in or¬ 
der to fully solve the analyticity conditions coming from 
the BES equations. Nevertheless their physical interpre¬ 
tation has not yet been addressed. The non-perturbative 
behaviour associated with singularities at s = nA (on 
the positive real axis), were seen to be directly linked to 
the mass gap of the 0(6) cr-model at least for n = 1 [28]. 
The singularities at s — —4 nA are much more suppressed 
and as such haven’t been addressed, but some questions 
can be immediately raised: does the same feature appear 
for the energy density of the 0(6) model? Does the rela¬ 
tion between scaling function and energy density hold for 
higher exponentially suppressed contributions? Since the 
two types of non-perturbative phenomena are collinear, 
will we witness resonance? 

For the aims of the current paper we used a transseries 
ansatz for the cusp anomalous dimension which did not 
include the second type of non-perturbative phenomena, 
as we only studied the solution up to n = 2 in the 
non-perturbative order. Nevertheless this was enough 
to check the resurgent properties of the strong coupling 
expansion, with the large-order relations predicted by 
resurgence accurately solving the large order behaviour 
of the perturbative expansion and first non-perturbative 
order. 

In our case we could determine the asymptotic expan¬ 
sions around the non-perturbative sectors via the BES 
equation, and use these to check the resurgence of the 
transseries. In cases when only perturbation theory is 
known, one can go the opposite way and use the predic¬ 
tions of resurgence to determine the coefficients of the 
expansions around the non-perturbative sectors. 

Knowing that the transeries proposed is indeed resur¬ 
gent we then proceeded to determine the resummed 
transseries, using a lateral resummation procedure. This 
naturally introduces an imaginary ambiguity, which can 
then be seen to cancel given the proper choice of the 
transseries parameter: using the methods of median re¬ 
summation [68] we determined the transseries parameter 
to be the one proposed in [28]. 

Finally the resummation procedure can be done for 


different values of the coupling, and we showed that in¬ 
cluding up to second order non-perturbative effects, we 
could systematically obtain accurate results for the cusp 
anomalous dimension all the way from strong coupling 
up to g ~ 0.2. Moreover we can perform the resumma¬ 
tion for any values of the coupling, as long as we take into 
consideration the Stokes phenomena occurring when we 
analytically continue our results across Stokes lines. 

This work is an example in the use and elegance of the 
resurgence techniques. Knowing the perturbative asymp¬ 
totic expansion of an observable in some regime, we can 
determine the non-perturbative corrections via large or¬ 
der relations, and upgrade the solution to a transseries. 
We can then use resummation methods such as Ecalle- 
Borel-Pade to obtain a resummed solution for any value 
of the coupling (using resurgence techniques to analyt¬ 
ically continue the solution across Stokes lines). This 
resummed transseries encodes the analytic properties of 
our observable as a function of the coupling. 

Other open questions still to be addressed are whether 
one can use a transseries ansatz for the auxiliary function 
in the BES equation, and through this find an equivalent 
approach to determining the coefficients r[ m ' ) to the one 
used in [24, 28]. Also, it remains to be understood how 
these results translate to the problem of the scaling func¬ 
tion e(g,j) and the energy density of the 0(6) cr-model. 
The first question is if one can find a full transseries in 
these cases. If so, will we have the same type of singular¬ 
ities? Finally, in these two cases we have two parameters 
g,j and it would be important to understand how the 
different regimes dictated by these parameters appear in 
the resurgent context. 


Acknowledgments 


I would like to thank Ricardo Schiappa, Romuald Janik 
and Hesam Soltanpanahi for useful comments and read¬ 
ing of an earlier version of this manuscript. This work 
was supported by the NCN grant 2012/06/A/ST2/00396. 


Appendix A: Asymptotic Expansions 


In order to write down the asymptotic expressions ap¬ 
pearing in the main text, we need first to define the fol¬ 
lowing four asymptotic expansions 


A 0 ,o (x,a) 
A 0 ,i {x,a) 


—'o r H) s! N s 

Y^ r Q + g)r(f+ s) 

ho r (!) s! (xa) s 


(Al) 

(A2) 
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Apo (x,a) 
Apr (x,a) 


v r(l + s ) 2 

ho r (!) s! (*«)*’ 

^ r(i + s ) 2 

ho r (s) s! M 8 ' 


(A3) 

(A4) 


In all the expressions in this appendix we assume a ^ 0. 
The functions [/„ (a;), n = 0,1 have been defined in [28] 
in terms of Whittaker functions of the second kind. For 
the purposes of this paper however, we are only interested 
in their asymptotic expansions for large x, which are: 


U o (y) ~ ( nx ) 5/4 A),i {x, -n ), 
u o (y) ~ ( nx )~ 3/4 A o,o (x,-n), (A5) 

U i (y) “ ^(nx) _ 1 / 4 Api(x,-n), 

U i (y) ~ \( nx )~ 3/i A hoi x ’~ n )- 


Note that even from these asymptotic expansions it is not 
difficult to see that these functions have a cut on the neg¬ 
ative real axis. The other functions of interest appearing 
in the main text are V n ( x ), n = 0,1. Again, in [28] these 
are entire functions which were written in terms of Whit¬ 
taker functions of the first kind, but for our purposes we 
are only interested in their asymptotic expansions. More 
specifically we will only need the asymptotic expansion 
of their ratio: 


2 r (a) = 2 



(A 6 ) 


The asymptotic expansion of this ratio for large x de¬ 
pends on the sign of a. For a = |a| > 0 


+oo 

2r (M) ~ 53(-l) fc e- fc l“l I e 5wifc / 4 (|a| x )~ k/2 R k (x, \a\) 


k—0 

(A7) 

where 


R 0 (x, |a|) 
R k (x, |a|) 


Apo (x, |a|) 
A 0 ,o (x, |a|) ’ 


(fIo (x, |a|) + |a| x 


Ayr (x, 

A 0 ,i (x, 



(AS) 

x 


f A 0 ,1 (x, — |a|) \ k 

V A o,o ( X , |a|) ) 


where 


R 0 (x, |a|) 
R k (x, |a|) 


Apr (x, |o|) 

A 0 ,i ( x , |a|) ’ 

Apo (x, 


Ro ( x , H) - 


A 0 ,0 (x, 


HA 

W\)J 


X 


(A10) 


( Aq,o (x, — |a|) \ k 
V A 0 ,i (x, |a|) ) 


In Section III when solving the analyticity conditions, 
some particular combinations of the asymptotic expan¬ 
sions Ai j (x, a) repeatedly appeared. These were 


(+)/ n _ A 0 ,i (x, —n) 


M^>(x) 


Ai i(x, —n) Apo (a;, a ) 

X I 1- 1 -V \OL\X+ -7-j—pr 

,A 0 i(a:,-n) A 0)0 (x, a ) 


+oo 


- H r i E x ~ r E 


M. 


(o,+) 


(All) 


r—0 


1=0 


M[ + \x) = 


A 0i0 (x, -n) 
xn (n + |a|) 


Apo (a:, -n) 
A 0 ,o (x, -n) 




-r-1 


M, 


App (x, |a|) 
Ao.o ( x , |a|) 
(i.+) 


r,£ 


r =0 


/ p | | r—i ’ 

e=o n \ a \ 


(A12) 


(-)/ N _ Aq.i ( X ,-n) 


M o ( x ) = 


(n+ |a|) 


a\ x 


x ( Api(ap -n) + Api (a:, joj) 
A 0 ,i (x, —n) A 0 ,i (x, |a|) 


+oo 


Mi Yt 


~r-i 


M. 


(o,-) 


r =0 


/ j o | | r—t ’ 

e=o n \ a \ 


(A13) 


m[ ^ (x) = 


A 0 ,o (x, -n ) 


a x 


x - 


Ap 0 (:r, -n) Apr (x, |a|) 


xn- 


A 0 ,o (x, —n) A 0 ,i (x, |a|) 


+oo 


- H r 7 5 > 


m: 


a-) 




r—0 


^=0 


n e |a| r 1 


(A14) 


For a = — | a | < 0 we then have 

+oo 

2r (- M) ~ J2(-l) k e- klalx e 3 ” ik/4 (|«| x ) k/2 R k (x, |a|) 

k—0 

(A9) 


Appendix B: Relations Between Sums 

Given the two fundamental objects ( n ) indepen¬ 

dent of the coupling x = 87 rg, appearing in the transseries 
of the coefficients c± (n, x) in (24), these can easily be de- 
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termined from analyticity conditions to be 


C +) (n) = 2 


r(n-f+ 1 ) 


(i) J 

^°’-)(n) = r(n ~^ + 2 l) . 

2n!r(f) 2 


(Bl) 

(B2) 


With these definitions we calcuate the following sums for 
to > 1 



3 

= g m+2'f'm+l ( a j l) j (B4) 


where p +\F p is a generalized hypergeometric function, 
a* are vectors with m + 2 entries: a + = (1, • • • ,1,5/4) 
and ad = (l,--* ,1,7/4), whereas b = (2, • • • , 2) is a 
vector with m+ 1 entries. Defining also A"q 0,± ' ) = 1/2, 
one can easily see that 


In all the expressions in this Appendix we assume ae = 
£—1/4, for f€Z. Other useful identities of 4°' (n) are 


E 

n> 1 


4°’ \ n - 1 11 ) 

n - \at\ 


= 0 , £< - 1 , 


(B6) 


i)£ 


n> 1 


4 ’ (n- 1) 

(n — \aie\) n m 


(B7) 


1 




K 


(o,-) 


^ M ro - s i«oi 


s— r+1 


, m,—£> 1, 


r l ? 


£ 

ri> 1 


4°’ (n-l) 


m ~ 1 


n" 


^=7, ™> 1, (B8) 

.£=0 l a °l 



4°’ V-1) 

(n — \ao\)n m 


m 


E 

t=0 


K 


(i,-) 


|«o| 


m—^ ’ 


TO > 0, 


(B9) 

where A'q 1, ^ = 7r/(4\/2) and for £ > 1 A"] 1- = 

(—1) X^l=o 4^ _a A'a°’ \ Useful identities of 4°’ + ( n ) are 


4°’ +) (n) |Q-| 
\4 J “-y (n — \a\)n m 

4o'~\n) |q-| 
V 4 / (n — \a\)n m 
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t=0 
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t=0 


K, 


(o,+) 


%-t ’ 


K, 


(o,-) 


I m—t ' 


(B5) 


b I V ^° ,+)(n) 

° n>i ( n + l a o|K 


= - E (-4)”*"* K, 


m—t jv-(l,+) 


1 =0 


where A', 


(b+) 


= 2 - 773 and K ‘ 


(i.+) 


= at; 


( 0 ,+) 


(BIO) 
for £ > 1. 
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